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We describe the bottomonium spectrum obtained on the UKQCD dynamical ensembles and its comparison to 
quenched results. We include a determination of a s and nib from the dynamical results. 



1. Introduction 

The T spectrum has long been recognised as a 
place to search for effects from including dynam- 
ical quarks because very precise calculations can 
be done and several mass splittings are sensitive 
to the short distance physics expected to be af- 
fected even by relatively heavy dynamical quarks. 
The calculations also enable the extraction of two 
important parameters of QCD, a s and mj,. There 
have been a number of calculations of the bot- 
tomonium spectrum using a range of dynamical 
configurations The results have not been 

conclusive, however. 

Here we give new results using dynamical con- 
figurations from the UKQCD collaboration [Q. 
These were obtained with Wilson (plaquette) glue 
and two flavours of clover dynamical quarks using 
a value of c sw given by the Alpha collaboration. 
Three ensembles are available with different dy- 
namical quark masses. These ensembles also have 
the advantage over previous calculations (for in- 
vestigations of this kind) of having the same lat- 
tice spacing. We believe that this enables system- 
atic effects as a function of the physical dynam- 
ical quark mass to be isolated more easily and 
not confused with effects that arise as the lattice 
spacing changes when the quark mass is altered 
at fixed beta. The dynamical quark mass varies 
from around m s to around 2m s and roughly 100 
independent configurations are available in each 
ensemble (separated by 40 trajectories). 

We use the standard 0(vf) NRQCD action for 



the b quarks || but we also include discretisation 
corrections to the E and B fields through the use 
of F^ v . This is given (tadpole- improved) by 
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The use of improved fields should reduce discreti- 
sation errors in the fine structure on these rela- 
tively coarse lattices ||. Two values of the bare 
quark mass rriQa were taken on each ensemble. 
They were chosen to give kinetic masses for the 
T particle which bracket the physical T mass, 
with the lattice spacing fixed from the splitting 
between the l 1 Pi and l 3 ^ states. This means 
that we can interpolate accurately to the physical 
T on each ensemble, again an important point in 
obtaining accurate results for the fine structure. 

We compare results on the dynamical ensem- 
bles to those on a reference quenched set at j3 = 
6.0, a slightly finer lattice spacing. 

2. Results for the spectrum 

Figure 1 plots the ratio of a -1 obtained from 
the IP-IS splitting in the T system to that de- 
fined from ro both on the dynamical configura- 
tions and the quenched set Q]. We have previ- 
ously checked that the ratio on quenched config- 
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Figure 1. The ratio of the inverse lattice 
spacing obtained from the IP-IS splitting of the 
T to that obtained from ro on quenched and 
unquenched configurations. The results on un- 
quenched configurations are plotted against the 
dynamical quark mass. 



Figure 2. The hyperfme splitting in MeV plot- 
ted for quenched configurations and (against dy- 
namical quark mass) for the dynamical configura- 
tions. Two values of heavy quark mass are given 
on each configuration and the physical b quark 
lies between them. 



urations is independent of the lattice spacing to 
within the statistical error. The results from the 
dynamical configurations are plotted against the 
square of the pseudoscalar light-light meson mass 
where the valence quark mass is equal to the sea 
quark mass. There is a possible sign that the dis- 
crepancy between the lattice spacing values ob- 
tained in the quenched approximation is reduced 
on the dynamical configurations. This is some- 
thing that we would hope for, because in the real 
world there is only one value of the lattice spac- 
ing. 

Figure 2 shows the hyperfme splitting between 
the T and the (so far unseen) r)b . In perturbation 
theory in a potential model this splitting is pro- 
portional to |^(0)| 2 , and is therefore dominated 
by very short distances. We then expect it to 
increase on unquenching as the screening of a s 
by dynamical quarks causes the potential at the 
origin to be deeper. This is clearly seen in our 



results, especially once the dynamical results are 
extrapolated to lighter values of the dynamical 
quark mass. 

The fine structure in the P wave sector is much 
noisier and no clear effect of unquenching can be 
seen. This is evident in Figure 3. 

3. Results for a s and m& 

We use the standard method of determining 
a s from a perturbative expansion of the plaque- 
tte fHU), setting the scale using the T IP-IS 
splitting. This requires a calculation of the Nf de- 
pendent part of the plaquette at 0(a 2 s ) for clover 
fcrmions. In the perturbative expansion to this 
order c sw — 1. We have (using also || and set- 
ting the dynamical quark mass to zero) 

1 4-7T 

- log < j-jrTrUp > = —a P {l - ba P ) 

a P = a P (3A0/a) (2) 
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Figure 3. Splittings among T P-wave states 
on quenched and unquenched configurations. 
Dashed lines give experimental results. 



Figure 4. The strong coupling constant deter- 
mined on various dynamical configurations at var- 
ious dynamical masses compared to that obtained 
on quenched configurations. 



with 

b = b 1 N c + b 2 N f (3) 

h = 0.39687 (4) 
121og3.40 - 10 
3d7t 

- 4ttP 4 + 32irx 2 (6) 

Pi = 0.006696001 - 0.0050467c™ (7) 

+ 0.029843(4, (8) 

x 2 = 0.00069292 - 0.0000202c sto (9) 

+ 0.00059624(4, (10) 



Figure 4 shows our results with the a s values 
run to a common scale of 8.2 GeV. Unqucnching 
shows up very clearly as o^ 2 ' is larger than . 
Comparison to other results shows a possible de- 
pendence on the type of dynamical quarks which 
we are investigating further. 

The b quark mass can be extracted using two 
different methods in NRQCD, either from the lat- 
tice bare quark mass or from the binding energy 
for, say, the T. To convert to the mass in the 
M S scheme we need either the lattice mass renor- 



malisation or the energy shift. Both of these are 
known only to 0(a s ) (TTJ at present so our re- 
sults here are not very precise. We see no effect 
from unquenching, nor do we see any variation 
with dynamical quark mass, see Figure 5. The 
most accurate results for mt, come at present from 
the binding energy of the B meson in the static 
limit 0. 

4. Conclusions 

Effects of unquenching are seen on the T hyper- 
fine splitting using the UKQCD dynamical con- 
figurations. The fact that they are matched in 
lattice spacing has helped with this. 

a s can be extracted using the experimental T 
spectrum and a perturbative expansion for the 
plaquette. This will give a value for a s to be used, 
for example, in perturbative matching of matrix 
elements to continuum QCD on the dynamical 
ensembles. 
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Figure 5. The b quark mass in the MS scheme at 
its own scale, determined from the binding energy 
of the T. Lattice and continuum perturbation 
theory to 0(a s ) are used. The results on the 
dynamical configurations are plotted against the 
dynamical quark mass. Two values are given on 
each ensemble corresponding to the bare lattice 
quark masses used. 
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